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Eq.  (2c)  should  be 

sa  V  =  T  -  0.7  Cn  PITS  -  eg  sin  9 

E2  U  O 

O 

line  16  of  text  should  be 

"The  values  of  Xt  and  X^.  obtained  .  .  . " 

line  24  should  be 

"is  used  as  the  estimated  value  of  Xt .  Tne  ... 
Eq.  (6)  should  be 
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wind  jgfects 

%  *  7=-*'71, 

V  sin  (T  -A) 

tan  £  *  - 

7n  +  v*  cos  (**  **A)  ' 


(5) 

(6) 


♦See  list  of  symbols  at  end  of  paper 

The  assumptions  underlying  these  equations  ai-e  veil  knovn  and  cay  be 
found  in  any  standard  text  on  airplane  stability  (e.g..  Reference  1,  Chap.  10), 

One  tacit  assumption,  however,  is  that  the  missile  is  roll  controlled  viih  no 
cross  coupling  between  directional  and  roll  maneuvers. 

The  simultaneous  solution  of  Equations  (1)  through  (6)  would  yield  all 
the  information  one  night  desire  about  the  flight  characteristics  of  the  nissile. 
However,  due  to  the  non-linear  nature  of  the  aerodynamic  coefficients,  analytical 
solution  is  difficult  if  not  inpossible  ana  one  is  farced  to  undertake  numerical 
or  analog  procedures  for  solutions.  The  rest  of  this  paper  Kill  outline  two 
approaches  to  the  complete  solution  of  the  equations  and  two  approaches  which 
by  caking  several  assumptions  reduce  the  equations  to  special  cases  for  specialized 
results.  Only  numerical  procedures  will  be  considered  since  analogs  involve  spe¬ 
cial  techniques  too  involved  for  this  limited  paper. 


A.  “Exact”  Analysis 

It  is  admitted  that-  any  solution  other  than  analytical  will  not  be 
exact.  However,  by  careful  and  refined  numerical  procedures  one  may  attain 
answers  which  my  be  as  close  to  the  exact  as  may  be  desired  for  engineering 
purposes.  In  view  of  this,  then,  exact  is  us id  here  to  indicate  that 
numerical  solutions  will  be  obtained  from  Equations  (1)  through  (6)  without 
making  additional  assumptions. 

There  are  several  numerical  methods  available  (Reference  2)  by  which 
one  may  obtain  solutions  to  differential  equations.  The  simplest,  ar.l  the 
one  most  often  used  in  boost-phase  analysis,  is  the  trapezoidal  rule,  the 
details  of  which  are  given  in  Appendix  I. 

Actually,  the  exact  analysis  is  seldom  undertaken  since  experience 
has  shown  that  two  simple  assumptions  lead  to  a  much  easier  and  faster  ap¬ 
proach  without  significant  loss  of  accuracy.  The  new  approach  will  be  dis¬ 
cussed  in  the  next  section. 


Semi-Exact  Analy 


The  semi-exact  analysis  is  very  much  like  the  exact.  However,  based 
on  experience  and  comparative  calculations,  it  is  assumed  that  oscillatory 
deviations  from  the  flight  path  do  not  affect  the  velocity  of  the  center  of 
gravity  along  the  flight  path  and  that  £  approaches  a  negligible  value 
quickly.  These  assumptions  allow  one  to  break  the  solution  of  the  equa¬ 
tions  into  two  parts:  the  determination  of  the  speed-time  history  alone 


"Ml  MUSS  UMtl'HJ 

waoswocmasn 
am  not  mtua 


Page  3 


and  the  solution  of  the  pitch  and  normal  force  equations  using  the  speed- 
tine  history.  This  particular  technique  was  used  by  the  Germans  in  their  V-2 
program  and  has  been  verified  numerous  tires  since. 


Assuring  that  oscillatory  changes  about  the  flight  path  do  not  affect 
the  forward  velocity  we  nay  rewrite  equation  (2)  (neglecting  the  wind)  as 


a  7n  *  T-D-cg  sin  '■C 

(2a) 

and 

since  at  t  *  D,  0  then 

^  90 

(ha) 

and 

b  Vn  "  (T"3)  -  eg  sin  90 

(2  b) 

The  procedure  is  tc  solve  equation  (2b)  by  the  trapezoidal  rule  to 
obtain  a  speed -tine  history  then  substitute  the  values  of  V  thus  obtained 
in  Equations  (l)  and  (3)  which,  taken  with  (It),  will  yield  a  complete  solu¬ 
tion.  The  results  are  as  accurate  as  the  exact  method  provided  care  is 
exercised  to  the  sane  degree. 

C.  KcDonneli  Speed-Tine  History 

Prcra  the  foregoing  section  it  must  be  clear  that  an  even  faster  or 
simpler  method  of  obtaining  the  speed-time  history  would  be  desirable. 
KcDonneli  Aircraft  Company  has  derived  such  a  method  in  Reference  3  by  making 
the  additional  assumption  that  speed  is  linear  with  time.  A  small  digression 
is  in  order  so  that  the  method  may  be  made  a  little  more  lucid. 

Equation  (2b)  may  be  written  fully  as: 

m  Vm  *  T  -  0.7  Cp  pM^S  -  mg  sin  0O  (2c) 

o 

where  the  constant  0.7  is  one-half  the  ratio  of  specific  heats.  Since,  in 
this  expression,  CpQ  is  a  complicat'd  function  of  Mach  number  a  correct  ap¬ 
proach  requires  the  trial  and  error  procedure  of  Appendix  I.  But.  recognizing 
that  the  speed-time  history  is  almost  lineal’  and  assuming  that  it  is  linear, 
one  may  compute  the  drag  for  any  time.  The  solution  of  (2c)  is  then  reduced 
to  a  summation  process,  and  this  is  the  heart  of  the  McDonnell  approach. 

To  apply  the  McDonnell  method  one  estimates  a  separation  speed  and 
time.  Then  assuming  that  the  speed  varies  linearly  with  time,  one  determines 
the  altitude,  pressure  and  speed  of  sound  at  each  time.  These  are  combined 
to  obtain  the  drag-time  history  and  Equation  (2c)  is  completed.  Details  are 
contained  in  Appendix  II. 

The  only  difficulty  in  applying  the  McDonnell  approach  is  in  esti¬ 
mating  the  separation  speed  but  a  reliable  method  for  this  has  been  developed 
and  will  be  discussed  in  the  next  section. 
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A  rel iable  estimate  Tor  the  separation  speed  is  required  for 
Section  lil-C  and  is  quite  often  necessary  in  evaluating  the  performance 
of  boosters.  A  simple  iterative  method  has  been  derived  by  KcCalley 
(Reference  Is)  based  on  an  approach  by  Walker  and  Henke  (Reference  5). 

The  equation  derived  by  KcCalley  is: 


where  p  and  Cg  are  selected  for  the  end  point  or  separation  condition. 

The  solution  of  this  equation  is,  then,  merely  an  evaluation  of  the  constants 
and  iteration  of  V  to  convergence.  The  details  underlying  Equation  (5)  are 
contained  in  Appendix  XU. 

IV,  Computation  Methods 

There  are,  at  AFL,  three  methods  of  computation  available  for  solution 
of  these  equations,  namely,  desk  calculator,  IBM  and  REAC.  Each  has  its  own 
peculiar  advantages  and  disadvantages  which  will  not  be  dwelt  on  here.  Instead, 
a  brief  note  as  to  which  is  adaptable  to  the  analyses  will  be  considered. 

The  REAC,  an  analog  computer,  could  be  used  with  the  exact  analysis  if 
enough  auxiliary  equipment  ware  available  but  since  such  side  equipment  is  not 
normally  available  use  of  REAC  is  restricted  to  the  semi-exact.  That  is,  if 
speed  histories  are  known,  then  the  rest  of  the  trajectory  may  be  worked.  It 
should  be  noted  that  even  this  requires  a  considerable  quantity  of  extra  equip¬ 
ment.  Since  the  bulk  of  this  paper  is  concerned  with  numerical  procedures  we 
shall  let  REAC  lie  on  this  brief  note. 

With  regard  to  the  numerical  procedures  it  is  true  that  IEM  can  do 
exactly  the  same  work  that  can  be  done  by  desk  calculators.  This  statement 
does  not  imply  that  it  is  true  a  priori  that  IBM  is  always  the  best.  There  is 
no  general  rule  which  may  be  invoked  to  determine  when  IBM  is  preferable  to 
desk  calculation  except  that  when  there  is  a  large  quantity  of  calculation 
IEM  is  essential.  It  should  be  fairly  obvious  the  IBM  machinery  would  be  un¬ 
called  for  on  the  McCalley  analysis  since  it  is  so  easily  worked  by  hand.  Other¬ 
wise,  eash  case  must  be  decided  as  it  arises. 


Accuracy 


When  obtaining  numerical  solutions  to  differential  equations  there  is 
always  some  loss  of  accuracy  from  the  possible  analytical  solution.  This  is 
more  especially  true  when  using  the  trapezoidal  rule  which  is,  from  a  mathematical 
standpoint,  rather  crude.  However,  if  care  is  exercised  and  incremental  changes 
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kept  snail,  the  errors  are  negligible  from  an  engineering  standpoint.  Therefore, 
it  will  be  assured  that  the  solution  to  the  “exact"  analysis  in  exact  and  the 
other  analyses  compared  with  it. 

The  serd. -exact  analysis  yields  practically  tha  sane  result  as  the  exact 
and  consequently  the  error  is  negligible.  In  fact,  a  check  was  made  against  data 
measured  in  an  actual  flight  to  check  separation  speed  only  and  the  difference 
between  calculated  and  ceasured  end  speeds  was  minute.  The  details  of  this  check 
are  contained  in  Reference  (6)  for  the  interested  reader. 

An  error  analysis  is  very  difficult  to  cake  cn  the  McDonnell  approach 
to  speed-tine  histories,  Numerical  comparison  has  been  cade  however,  and  the 
results  show  that  the  difference  between  the  send. -exact  and  the  McDonnell  speed- 
tine  history  is  generally  less  than  2$  when  the  estimated  separation  speed  is 
between  2  and  h%  in  error  (Reference  7).  For  cost  trajectory  analyses  the 
McDonnell  speed-time  history  is  sufficiently  accurate  since  a  snail  error  in 
velocity  at  any  given  time  does  not  appreciably  affect  the  results. 

The  HcCalley  method  of  computing  separation  speeds  produces  errors 
up  to  h.%  depending  on  the  Mach  number.  It  appears  (Reference  7)  that  the  error 
increases  with  increasing  Kach  number,  which  is  not  surprising  since  the  method 
assumes  time-linear  drag.  Hence,  the  higher  the  separation  velocity  (Kach 
number)  in  a  given  time  the  greater  the  drag  deviates  from  a  straight  line  and 
the  greater  the  error. 

VI.  Application  * 

How  that  the  various  methods  are  in  hand  it  would  seem  that  a  statement 
on  the  uses  of  each  would  be  in  order.  Experience  has  indicated  the  area  of 
usefulness  to  which  each  is  suited  and  the  types  of  problems  which  have  been  en¬ 
countered. 

As  must  be  evident,  the  exact,  or  since  the  difference  is  small,  the 
semi-exact  analysis  is  called  for  whenever  the  interest  is  in  the  actual  flight 
performance  of  the  missile.  Either  of  these  will  yield  time  histories  of  all 
of  the  variables  in  question  which  are  essential  in  computing  air-loads,  inertia 
loads,  dispersions,  etc.  This  field  of  requirements  is  self-evident  and  usually 
no  question  will  arise  as  to  when  to  use  the  full  analysis. 

If  one  is  interested  only  in  the  speed-time  history  then  some  conflict 
appears  as  to  the  better  method.  In  fact,  even  when  making  the  semi-exact  ana¬ 
lysis,  the  question  arises  as  to  whether  an  approximate  speed-time  history  such 
as  obtained  from  the  McDonnell  method  is  sufficient  for  trajectory  analysis.  It 
is  the  writer’s  opinion  that  it  is.  This  opinion  is  based  on  the  fact  that  the 
aerodynamic  coefficients  do  not,  throughout  the  largest  portion  of  the  boost- 
phase,  change  radically  with  Mach  number.  Thus  if  cne  has  a  small  speed  error  at 
some  step,  it  is  obvious  that  the  forces  computed  are  not  greatly  different  from 
the  exact.  The  error  generated  thusly  is  more  than  compensated  for  by  the  saving 
of  time  and  expense.  However,  one  must  judge  for  himself  on  this  point. 

There  are  times  when  the  separation  speed  is  of  the  utmost  importance. 
For  example,  a  certain  minimum  speed  is  required  for  engine  starting  and  there  may 
be  a  question  as  to  whether  the  booster  impulse  is  sufficient  to  produce  this  s^eed. 
In  such  a  problem  the  semi -exact  speed  history  is  an  absolute  must. 
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Whenever  the  problem  of  booster  design  arises,  there  are  several  ways 
to  attack  it.  First,  by  using  the  McCalley  method  one  may  determine  a  general 
area  of  impulse  necessary  to  reach  a  given  speed.  The  impulse  may  then  be  trans¬ 
lated  into  more  exact  speeds  by  using  one  of  the  more  refined  methods.  Further¬ 
more,  one  may  show  the  effects  of  changes  in  the  booster  and  missile  parameters 
by  using  the  McCalley  equation.  Thus,  in  an  initial  design  study  the  KcCalley 
approach  is  extremely  useful  in  setting  approximate  boundaries.  All  of  this 
presupposes  that  the  propellant  characteristics  are  known  or  estimable  as  they 
relate  to  physical  dimensions,  gas  pressures,  etc.  If  one  is  starting  from 
scratch  then  a  more  refined  approach  such  as  that  developed  by  Hawley  and  Fenton 
in  Reference  8  is  required. 

VH.  Techniques 


There  are  certain  details  which  arise  during  the  actual  calculations 
which  are  sometimes  confusing.  This  section  will  try  to  clarify  a  few  of  the 
more  important  oner ; 

A.  First  consideration  will  be  given  to  wind  effects  and  how  these  are  included 
in  the  various  analyses.  When  using  the  exact  analysis  one  must  account  for 
the  wind  at  every  step.  Obviously  this  amounts  to  solving  Equations  (1) 
through  (6)  without  modification.  It  should  be  apparent  that  the  relative 
velocity  %  is  used  in  computing  the  aerodynamic  terms.  To  illustrate,  take 
Equation  (2)  and  write  in  full:  ^ 

m  Vm  «  T  cos  (c /-  £)  -  Cp  S  cos  c  -  mg  sin'/’  (6) 

Thus,  the  Vra  refers  to  acceleration  relative  to  earth  while  the  refers 
to  missile  velocity  relative  to  the  air. 

The  semi-exact  method  may  fortunately  be  handled  a  little  differ¬ 
ently.  Since  the  speed  of  the  missile  becomes  much  greater  than  the  >dnd  in 
a  very  short  time  then  one  is  justified  in  neglection  of  the  angle  £.  . 
Therefore,  the  only  things  one  must  consider  is  the  change  in  angle  of  attack 
which  is  defined  by 


V-y  sin  (9  ~X) 
tan  <K -  Vm  +  Vw  cos  (9  -A) 


(7) 


where  positive  Vw  is  a  headwind  in  this  scalar  equation. 

In  the  McDonnell  analysis  one  must  assume  a  velocity  history  based 
on  an  average  acce1  sraticn.  To  compute  the  drag  the  relative  velocity  is 
used  so  that  one  must  add  at  every  time  the  longitudinal  component  of  the 
wind  to  the  speed-time  history  computed  from  the  average  acceleration. 

The  McCalley  equation  makes  no  provision  for  wind  but  it  may  be 
considered  that  the  equation  always  yields  the  relative  velocity  as  the 
solution.  That  is,  the  V  obtained  by  iteration  already  contains  a  wind 
component. 
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B.  A  technique  often  used  by  McDonnell  is  that  of  calculating  everything  in  the 
so-called  "slant  plane".  This  amounts  to  rotating  the  coordinate  axes 
through  the  angle  90  or  the  original  launch  angle.  There  is  some  simplifica¬ 
tion  achieved  in  this  way  and  in  fact  by  neglecting  gravity  one  may  determine 
a  sort  of  "universal"  trajectory.  The  principal  advantage  is  in  the  relative 
size  of  the  numbers  handled.  The  details  are  not  important  since  all  that  is 
required  is  rotation  of  the  axes. 


C.  In  setting  up  the  McDonnell  analysis  for  solution  on  IBM  equipment,  one  en¬ 
counters  the  difficulty  of  non-linear  pressure-altitude  or  density-altitude 
relations.  If,  however,  these  relationships  are  expanded  in  a  Taylor  series 
about  the  altitude  of  the  launching  site  then  by  taking  the  first  two  terms 
the  pressure  (or  density)  nay  be  determined  approximately  in  linear  form. 

It  has  been  found  that  no  substantial  error  is  introduced  up  to  an  altitude  of 
about  3>000  ft  above  the  launching  site.  The  same  sort  of  approximation  may 
be  made  for  the  speed  of  sound. 


D. 


In  working  the  speed  histories  using  the  semi-exact  analysis  and  the  trape¬ 
zoidal  rule  it  has  been  found  that  a  sizeable  error  jb ay  be  allowed  in  the 
estimated  value  of  V  (or  as  notated  in  Appendix  I,  X).  The  reason  for  this 
is  apparent  through  examination  of  the  following: 

9  • 

Let  X-t  be  the  correct  value  of  the  acceleration  at  the  end  of  the 
next  time  increment  and  Xg  be  the  estimated  value.  Then 


*t 


•  •  •• 

xt  +  At-1 
2 

If  •# 

xt  +  XE  *  R 
2 

•  •  #f 

+  xE  .  .  R  A  , 

- x  —  A.  t  +  It.!  +  2  At 


At  +  Xt-1 


A  t  +  Xt.! 


Xt 


9  m 


Xt  +Xt-1 


At  +  x. 


t-i 
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-A-— ■  -  A  t  +  2  Xt-1  +  2  A  t 


A  t  +  x 


t-1 
R 


xt  *  Xg  "  o 

— : -  A  t2  +  Xt.!  A  t  +  Xt„!  +  ;  Ar 


„  This  shows  that  the  value  of  Xt  is  in  error  by  the  amount  of 

^  At2.  If  At  =  0.1  sec  and  R  *=  100  f£/sec 2  then  the  error  in  Xt  is 

only  0.2£  ft  with  an  even  smaller  change  in  altitude.  The  error  in  the 
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speed  L.  is  2  At  ?cd.ch,  nsirg  tbs  ssre  rzabers  as  before,  is  5  ft/sec. 
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?fersrical  Integration  by  Trapezoid  Rale 


Assure  the  following  equation  is  to  be  solved 

X  -  b£  +  el  +  d  (1) 

where  the  fe,  e,  and  d  are,  in  general,  functions  of  X  and  t,  and  that  we  have 
reached  the  solution  of  the  equation  up  to  sore  tire  t  *  '£. 

5y  extrapolating  the  past  history,  or  by  any  Keans  available,  one  esti¬ 
mates  the  value  of  X  after  sore  small  tire  increrent,  (t  -  T),  as  X^.  If  the  tine 
increment  is  srall  enough  and  the  functions  behaving  sroothly  then  an  average  value 
of  X  nay  be  calculated  for  the  interval  as: 


It  +  Itr 


avg 


Then  the  value  of  X  *  X^  at  the  end  of  the  interval  is  approximately 


It  “  lavg(t-t)  +x^ 


•  #  #• 

it  +  i-c 


(t  -  T)  *  ZZ 


and  the  value  of  X  »  Xj.  is 


It  “  x>vg  (t  -t)  +  X- 


>vg 

L  +xT 


(t  -T)  +  X-c 


(2) 


(3) 


(U) 


The  values  of  X^  and  Xfc  obtained  from  Equations  (3)  and  (Ii)  are  then  multiplied  by 
b  and  e,  which  are  selected  for  t  *  t  and  x  *  xt  and  added  to  the  value  of  d  # . 

„ elected  the  same  way.  The  sum  thus  obtained  called  the  calculated  value  of  X^ 
is  compared  to  the  estimated  value  of  X-^  and  if  agreement  is  attained  one  proceeds 
to  the  next  time  increment.  If  agreement  is  not  realized  then  X-^  is  estimated  again 
until  the  estimated  and  calculated  values  of  are  very  nearly  the  same. 

If  in  Equation  (l)  b,  e  and  d  are  constants  or  are  functions  of  t  only 
then  the  process  becomes  iterative,  i.e,,  the  value  of  X+  obtained  from  calculation 
is  used  a3  the  estimated  value  of  Xt.  The  process  is  usually  convergent.  Obviously, 
if  b,  e  and  d  are  constants  then  an  analytical  solution  is  straightforward. 
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APPENDIX  II 


McDonnell  Method  of  Computing  Speed-Time  Histories 


as: 


Assume  the  following: 

1.  The  speed-time  history  is  linear 

2.  The  separation  speed,  7S 

3.  The  separation  time,  ts 

Hith  these  assumptions  one  proceeds  by  calculating  an  average  acceleration 


aavg 


(1) 


Then  compute  the  speed 

V(t)  =  aaVg  t  +  V0  (2) 

and  the  slant  range 

S(t)  -  2  aaVg  t^  +  VQt  +  SQ 

The  altitude  h,  is  then  calculated  from  S(t)  and  the  speed  of  sound  and  pressure 
determined  as  functions  of  h  as  follows: 


h  *  S(t)  sin  -0O 
c  -  fi(h) 

P  =  f2(h) 


(3) 

(U) 

(5) 


steps: 


The  drag,  D,  as  a  function  of  time  is  then  easily  found  by  the  following 

V(t) 

1.  Compute  the  Mach  number  M(t)  =  c 

2.  Select  the  drag  coefficient  corresponding  to  M(t)  at  each  time 
interval. 

3.  Complete  the  drag  calculation  by  the  following  equation: 


D(t) 


0.7  Cq  IPS 


(6) 


where  the  0.7  is  one-half  the  ratio  of  specific  heats  for  air. and  S  is  the  reference  area 


The  remainder  of  the  speed  history  is  just  a  question  of  evaluating  the 
following  equation  and  applying  the  trapezoid  rule. 


a(t) 


T(t) 

Ml'  + 


D(t) 

MTtT 


g  sin  9 


(7) 
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APPENDIX  IH 


s 

t  HcCalley  Equation 


The  KcCalley  equation  is  based  on  a  method  developed  by  Walker  and  Henke 
wherein  the  drag  is  assumed  to  be  linear  with  time.  Under  this  assumption  plus  the 
additional  one  that  the  mass  varies  linearly  with  time.  Walker  obtains  the  follow¬ 
ing  equation:  /  " 


G 

2 


G 


I 

2 


u 

(g) 


Tt, 


% 


-In 


(1 


+  gt  sin  -9 


"I 


(1) 


where 


(2) 


McCalley's  contribution  comes  about  by  restating  the  Walker  equation  as 


(3) 


and  then  expanding  the  logarithm  terms  in  a  series.  Taking  only  the  first  two 
(linear)  terms  and  applying  correction  factors  McCalley  derives  the  following  ex¬ 
pression: 


T  t. 


%  “  2  % 


Kt  -  g  t  sin  6 


*n+l 


s  t. 


U  (Mo  -  3  %) 


K2  Vn 


Ot) 


Observation  of  drag  calculations  made  under  more  exact  conditions  has  led 
to  the  conclusion  that  the  drag  is  linear  not  from  t  ■  0  but  from  t  -  1  sec.  This 
fact  in  turn  leads  to  a  modification  of  Equation  (U)  which  is  more  nearly  in  keeping 
with  the  physical  situation: 
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SYMBOLS 

Aerodyr  amic  Coefficient 
Acceleration,  ft/sec^ 

Aerodynamic  Coefficient 
Aerodynamic  Coefficient 
Aerodynamic  Coefficient 
Speed  of  sound,  ft/sec 
Drag  force,  lbs 
Aerodynamic  Coefficient 
Aerodynamic  Coefficient 
General  Coefficient 
Aerodynamic  Coefficient 
General  Coefficient 

Defined  in  Equation  (2)  of  Appendix  III 
Gravitational  Constant,  ft/sec^ 

Altitude,  ft 

f-^  (Mp/M0),  Correction  Constant 
f2  (Mp/M0),  Correction  Constant 

Perpendicular  distance  from  center  of  gravity  of  missile  to  action 
line  of  thrust,  ft 
Mach  Number  ■  V/C 
Mass,  slugs 

Total  Mass  at  Lai'nch,  slugs 
Mass  of  useable  propellant,  slugs 
Atmospheric  pressure,  lb/ft^ 

Difference  between  assumed  and  true  acceleration,  ft/sec^ 

p 

Reference  Area,  ft 
Slant  Range,  ft 
Thrust  Force,  lbs 


t  i 
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SIMBOLS  (Continued) 


Page  2 


Velocity,  ft/sec 

Distance  measured  a3  . ng  launching  line,  ft 
Zero  angle  of  attack  drag  coefficient 
Angle  of  Attack 
Flight  Path  Angle 

Surface  deflection  either  by  control,  deformation  or  manufacturing 
difficulties 

Angle  between  relative  wind  and  flight  path 
Angle  between  missile  centerline  and  reference  plane 
Angle  between  reference  plane  and  wind  vector 
Air  Density,  slugs/ft^ 


Subscripts 


Estimated  Condition 
Iteration  Number 
defers  to  Properties  of  Missile 
Initial  Conditions 

Refers  to  Conditions  relative  to  Missile 
Separation  Conditions 
Conditions  at  Time  t 

Conditions  at  Preceding  Time  Considered  (not  necessarily  a  unit 
increment  of  time  such  as  1  sec) 


auperscnt 


Dots  -  Refer  to  derivatives  with  respect  to  time  e.g,,  X  =  ” 

dt 


Bars  -  Refer  to  vector  quantities 
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